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AIRFOIL THEORY AT SUPERSONIC SPEED* 
By H. Schlichting 



A theory is developed for the airfoil of finite span 
at supersonic speed analogous to the Prandtl airfoil theory 
of 1918-19 for incompressible flow. In addition to the 
profile and induced drags, account must he taken at super- 
sonic flow of still another drag, namely, the wave drag, 
which is independent of the wing aspect ratio. Both wave 
and induced drags are proportional to the square of the 
lift and depend on the Mach number, that is, the ratio of 
the flight to sound speed. In general, in the case of 
supersonic flow, the drag-lift ratio is considerably less 
favorable than is the case for incompressible flow. Among 
others, the following examples are considered: 

1. Lifting line with constant lift distribution 

(horseshoe vortex). 

2. Computation of wave and induced drag and the twist 

of a trapezoidal wing of constant lift density, 

3. Computation of the lift distribution and drag of 

an untwisted rectangular wing. 

£ 

I. INTRODUCTION 



The basic principles for the following computation of 
airfoil flow at supersonic speed are presented in the paper 
of professor Prandtl (reference 1) , and a detailed expla- 
nation of the method may therefore be dispensed with here. 

The potential $ ^ of a lifting line at supersonic 

speed may be derived in a simple manner from the potential 
$q of a stationary source in the presence of a supersonic 

flow. 

* "Tragf lugeltheorie bei Uber schallgeschwindigkei t . 11 Jahr- 
buch 1937 der deutschen Luf tf ahrtf or schung, pp. I 181-97. 
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If $q denotes. the source potential of strength 4tt, 
the potential $ T of the lifting, line element dy with 

circulation T about the y axis is given by 

x' = x 

* T - a — /* * 5_fia (i) 

1 4 TT J d Z 

The potential of a source at the point x = y = z = 0 in 
the presence of a flow, with velocity u Q > c in the di- 
rection of the positive x axis is 

.«»■■■ 7 . 1 ===== = >. 

/x a - [(— ) 2 - x ] (y 2 + zS ) 

The potential (2) is real within the double cone with half 
cone angle a, the axis of which cone is parallel to the 
direction of flow (sin a= c/u 0 ). Outside of this cone 
the potential, according to the formula* is imaginary, 
• Actually, $q is there to be taken identically equal to 
zero. The potential has physical reality only in the "af- 
ter cone" of the point x = y = z = 0, In the "forward 
cone" it is similarly to be taken identically equal to 
zero, 

The potential $^ is the starting point for con- 
structing the airfoil potential. We shall first derive 
from it the potential of a line source of finite length, 
then with the aid of the operation indicated inequation 
(1) we shall obtain the potential of a lifting line of 
•finite length for various lift distributions. From the 
lifting line, there is finally obtained by the familiar 
method/ the lifting surface. In this manner,, a theory 
of the airfoil of finite span for supersonic speed is ob- 
tained that forms the counterpart of the Prandtl airfoil 
theory for the incompressible flow case (reference 2). 
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II • CONSTANT LIFT DISTRIBUTION (LIFTING- LINE) 



We shall now assume a line source of length h 
(later = span of wing) which lies in the direction of the 
y axis and extends from y 1 = -b/2 to y* = +b/2 (fig. 1). 
Let g(y f ) he the initially given local source intensity 
(later = the lif t di str ihut ion) . Further, let x f y t z 
he the coordinates of a point in the flow and 0, y ! , 0, 
the coordinates of a source point. Then from equation (2) 
the potential of the line source is 



y* = + l 



g(y') d y' 



(3) 



We introduce nondimensional coordinates hy dividing all 
lengths hy the half-length h/2 of the line source and 
accordingly set 

. t lt b ^ v & * b 71 

Further, we introduce the abbreviated notation 



or 




where a denotes the Mach angle. 

The potential of the line source then becomes 



r,«= +1 



g(*n') d <n» 



(5) 



^4 '9 ? K f " 
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We shall now carry out the Integration in equation 
(5) for the simplest case Of the line source with constant 
source density, that is t - g-CO*') = const. = 1. This gives 

d T}« 

»a = / ■ (6) 



Writing T| - T|» = o; T>» = -1; o a d x » n + 1 

T)» = +li 0 = * 2 = T) - 1 

and (I/*)* - I* = a 8 

equation (6) becomes 

d = 3. 



By the operation in equation (1) there is then obtained 
the potential $ T of the lifting line with the constant 

lift distribution r o , setting 



4 o / » ° *0 



(3) 



The first step of the above operation t differentiation with 
respect to i $ may he carried/out immediately but the inte- 
gration requires a somewhat longer computation. There is 
obtained 

- K —2 = 3 - (9) 

95 aVa 3 - * a 2 a* 7a* - V 

In integrating with respect to I, $ a and $ x are con- 
stant, For the first term, there is obtained 
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p V r £ * s *e* , (lQ) 

where there has "been set 

t* = e - u t) a 

The evaluation of the integral gives 



a 2 7a 2 - * 2 2 



- 2 k 2 * 2 £ 2 + U* - * 2 8 ) V 



- — arc tan 



2 2t« B /i* 2 + k 2 e*(t 2 - v> - K 4 £ 2 * 2 



'2 



£ 2 l s (-n - i) a 

= — arc tan — (11) 

2 2 ! (ti - 1) t-/uT" 

where 

ui = i s - k 2 [(n - i) 2 + t 3 ], (12) 

Since the integral (10) outside of the Mach cone, 
at the end point T\ = 1 of the lifting line with axis 
parallel to the x axis, 

I 2 - k 2 [(ti - I) 2 + t*l = 0 

is imaginary, i.e., is to he taken equal to zero, the in- 
tegration with respect to | f need not he extended from 
I 1 = -oo but only from the cbne surface along lines paral- 
lel to the x axis. For the lower limit of integration, 
we have thus the constant arc tan coiTT/2, which we may 
suppress. There' is thus found from (8), (9), and (LI) the 
required potential of the lifting line with constant cir- 
culation 
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i a u, - e 3 ^ - i) 2 



e 1 = t . 

- — * « / ' ? 11 = - — arc tan 

• 2tt J a^Va 8 - d 3 4tt 2 t '(T) - l)t </TiT 

e »=-» 



3 



r 0 I (T| - X) 

= — arc tan + similar term for the cone at 



This potential is different from zero only within the two 
iach cones arising at the ends of the lifting line (u) > 0) 
while in the entire remaining space it is equal to 0. For 
a complete circuit about each of the cone axes *n = ±1, 
i = 0 f the arc tan increases by 2 tt. The enclosed vortex 
filament therefore has the circulation P 0 . The lifting 
line assumed to extend from y = -b/2 to y = +b/2 with the 
constant circulation P 0 along the span continues behind as 
a free vortex .line in ".the two axes of the Mach cones* Equa- 
tion (13) thus gives the potential of a "horseshoe vortex* 
at supersonic flow. As in the case of the incompressible 
flow, this simple horseshoe vortex becomes the starting point 
for more complicated lifting systems* 

In order to obtain an idea as to the appearance of the 
supersonic flow in the neighborhood of a horse&oe vortex, 
we differentiate the potential (13) to find the induced ve- 
locities 



C v = - 



3$ 




. a z 



and obtain 
c„ 



= K 



C- = - 



On - i) t 



tt b (f - k t 8 ) vnr 



it b C i"+ - D 3 ] /uT 

eh - i) k*i 2 ) 

tt b 'it B „+ (n - i) a i 8 ] jz 



> 



(14a, D,c) 



x See footnote on next page. 
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The field of these velocities exhibits a number of singu- 
larities. On the cone surface all three velocity compo- 
nents become infinite. On the cone axis c x = 0, hut 
c y and c z become infinite as l/r (where r is the 
distance from the axis). In the neighborhood of the cone 
axis, c y and c z thus behave exactly as in the neigh- 
borhood of a vortex filament in the incompressible flow. 
The field of the induced velocities gives a motion which 
encircles the vortex filament traveling downstream from 
the end of the lifting line Tj =1, |a.t = 0, as may 
be seen immediately from (14). 

In the plane T| - 1 = 0 through the end of the lift- 
ing line c 2 = 0 and 

i > 0 : c y < 0 

£ < 0" : .c y > 0 

In the plane t = 0 t which contains the lifting line, 
Cy = 0 and 

Ti-l>0:c z >0 
T) - 1 < 0 : c 2 < 0 

The flow picture in the cone, however, in its* detail is 
essentially different from that in the neighborhood of a - 
vortex filament in the incompressible case. Figure 2 shows 
the flow picture of the y and z velocities in a plane 
perpendicular to the cone axis downstream of the lifting 
line. The figure was obtained by computing the isocline 
field c 2 /cy = const. 'On the cone surface, as has been 
said, c 2 and c y are infinite, although for the slope 
of the streamlines- -c 2 /c y there is here obtained the sim- 
ple value 



c 2 _ t 

Cy T| - 1 



*A check for the correctness of this solution is obtained 
by substituting in the linearized continuity equation 



- K a 



8 x 




which must be identically satisfied 
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The direction of the streamlines is therefore radial to 
the center*- The flow consists partly of the closed stream- 
lines which circulate about the vortex filament and partly 
of the streamlines that enter on. one side of the cone and 
leave it again on the other side. 

In. addition to. the two Mach cones that arise froA 
each of its ends, the lifting line generates two plane 
waves, which enclose a "wedge space" and which appear in 
the streamline picture as the common tangents of the two 
cones. 

For the downwash distribution in the plane i = 0 
through the cone center, there is obtained from (14c) the 
simple formula 

2 tt x J\ - d 2 

~ tan a c 2<> - (15) 

o 

where. 

b 

y 

4 = (15a) 

x tan a 

This downwash distribution is shown in figure 3. 

In order to study the processes on an airfoil of fi- 
nite length at supersonic speed, particularly the induced 
drag, the replacement of the wing by a lifting line with 
constant circulation as in the case of the incompressible 
flow, appears inadmissible since on account of the infi- 
nite velocity at the end of the lifting line an Infinite 
induced drag would be obtained. This difficulty in the 
case of the incompressible flow is avoided, as is known, 
by allowing the circulation to drop to zero in a suitable 
manner toward the wing tips. The induced drag is then 
computed by the formula. 

y = +b/e 

>I = P J o Zo (y) r (y) d y (16) 
y = -b/ 8 



(where c z Q (y) ls induced downwash velocity at the 

place of the lifting line, and P the density). 
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In the case of the supersonic, flow, the relations are 
complicated by the fact that in spite of the assumption of 
a lift distribution decreasing to zero toward the wing tips, 
there are obtained singularities at t he lifting line posi- 
tion of such a character as to make the computation of the 
induced drag by formula (16) f which maintains its validi.ty 
for supersonic flow, impossible. As closer investigation 
shows, this is due to the fact that the lifting line is the 
geometric locus of the verticfes of all the" Mach cones that 
pass down behind. This difficulty may be overcome by pass- 
ing from the lifting line to the lifting surface. 

III. WAVE RESISTANCE (DRAG ) 

Before proceeding to the corresponding computations, 
we shall discuss briefly the supersonic flow about an in- 
finitely long airfoil (two-dimensional problem), a problem 
that had been considered by J. Ackeret in 1925 (reference 
3). 

The simplest and at the same time the ideal supersonic 
profile is that of the infinitely thin flat plate of chord 
t set to a small angle of attack p Q (fig. 4), For such 
a plate the lift per unit span is 

A = 2 tan a $ 0 t P u. s (17) 



or 



A 



P 2 * 
r u rt t 



= c a = 4 tan a £ ft (18) 



2 



On account of A = p u Q T Q , the relation between the an- 
gle of attack of the wing and the circulation is 

1 T 

8 A u A tan a = (18a) 

P ° 0 2 t 

From the incompressible flow, the supersonic flow 
about the airfoil differs in that,- for the latter case, 
even if the fluid friction is neglected, there is always 
associated a drag that originates from the plane waves 
which start out from the lifting surface and are inclined 
to the latter by the Mach angle and which, therefore may 
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be denote* as the wave drag. For the flat plate, the wave 
drag per unit span; is 

*wave = - P 0 A = - 2 tan a Po 2 t P" V C 1 *) 

or 

— — = c w = 4 tan ap 0 3 (19a) 

r tt. t 
2 .0 

The resultant of the lift and the drag Is here at right 
angles to the plate. This comes from the fact that at su- 
personic flows there is no suction force at the leading 
edge of the plate. From equations (1$ and (19a), there is 
obtained for the polar of the wave drag 



c w 



c 2 
c a 



wav * 4 tan a 

which is thus a parabola as in the case of the incompres- 
sible flow. 

Plane waves start out from the leading and trailing 
edges of the inclined flat plate (fig/ 4) and in the space 
between them the induced downwash v elocity is 



"wav« ■ o » . 2 t tan a 



c Z __ = -P fl u 0 = - i 2 (20) 



The wave drag, on the other hand, can also be computed from 
this downwash velocity induced by the plane waves, accord- 
ing to the formula 

*wave = P r b c *wave (21) 

as may be seen by comparison with (19) and (20). In the 

next section it will be shown that, for a lifting surface, 

the velocity induced by the tip vortices like c„ is 

. z wave 

proportional to T 0 /t tan a. It then follows from equa- 
tions (21) and (16) that the wave drag behaves in exactly 
the- same way as the induced drag from thw. tip vorticos 0 .. 
For practical applications it is tharefore of no inter- 
est to consider the induced drag alone, but it is the sum 
of the induced and wave drags that must be considered. 
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For an airfoil of finite span and constant chord with 
circulation that is constant along the chord and variable 
along the span T(y) = t 7 (y) the total lift and wave 
drag are given by 

y = +t>/2 r,= +.i 

A=Pu o y r d y = i p u 0 b t J y (r\) d -n (22) 

y = _b /a T)= -i 

+b/2 r\ = +1 

W wave = p f c *ow T d 7 = § b * f c *o w 7 d ^ (23) 
y = -b/2 T| = — x 

where 

1 7 



° w 2 tan a 

is the induced wave velocity. Accordingly 



(24) 



"wave = " " / 7 2 d 11 (25) 

4 tan 



i a J 

•n = -i 

By comparison of equations (22) and (25), there is found 
the relation between drag and lift 



W - 9 1 * / A \2 

w wave - 2 - 



(3.) 

Nil -h/ 



p t tan a x u 0 b 

In the above equation Z is a nondimensional coef- 
ficient that depends only on the lift distribution 



(27) 




From equation (26), it follows that: 
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'wave 



tan a 



-a 



(28) 



The numerical values of Z are given in table I for sev- 
eral simple lift distributions. 



TABLE I - Values of the Coefficient Z for Various 
Lift Distributions (Wing of Rectangular Plan Form) 



Number 


Lift distribution 


Z 


1 


rectangular 


1/4 




= 0.250. 


2 


ellipse 


8 

3 TT 3 




= .271 




trapezoidal 
« b« = b/2 


2 


b 




3 






2 
) • 

= .370 


4 


parabola 


-3_ 

10 




= .300 


5 


triangle 


2 
3 




= .667 



IV. LIFTING SURFACE WITH CONSTANT LIFT DISTRIBUTION 



For the successful computa 
for supersonic flow f according 
taneous assumptions must be mad 
lift toward the edges and a sur 
bound vortices. This twofold e 
a considerable swelling of the 
of induced velocities as compar 
flow where the computation invo 
In order to be able to recogniz 
of each of these two extensions 



tion of the induced drag 
to section II, the simul- 
e of a suitable drop in 
face distribution of the 
xtension means naturally 
computation of the field 
ed with the incompressible 
Ives mostly a lifting line, 
e more clearly the effect 
we proceed in two steps. 



N.A.C.A. Technical Memorandum No. 897 13 



We first maintain the lift distribution constant along 
the span and consider only the transition from lifting 
line to lifting surface. The field of induced velocities 
thus obtained for a wing with constant spanwise circula- 
tion distribution and constant chord,: while it does not 
enable as yet the computation of the induced drag never- 
theless furnishes many useful results so that we proceed 
first to compute this field. 

We assume therefore the. circulation T 0 constant 
along the Span b as uniformly distributed over a rec- 
tangular lifting surface of chord t and extending from 
x = 0 to x = t (fig. 5). The circulation for a strip 
of the lifting surface of unit width is therefore 7 0 = 
^oA» I* would be most convenient to make the transition 
from the lifting line to the lifting surface directly on 
the potential (13), On account of the integration diffi- 
culties that arise, however, the transition will be made 
on the velocity components (equation (14)), first for the 
z component since the latter is the most important for 
the computati on of the induced drag. 

A strip of the lift ing surf ace of width dx 1 at a 
distance x* from the leading edge contributes to the 
induced z component . c z at the point x, y, z, if the 
point lies within the Mach cone arising from the end of 
the strip the amount 

d ? z = ^2- d x' f(£ - f.TiJ) = -I 0 - d £• f a - I'.Tl.O 

TT t) . 2 TT ; 

where, according to equation (14) 

f(M.D- ^-»<»-'t') (29) . 

U 2 u> + (r\ - l) 2 ? 3 ]^ 

If the point x, y, z lies outside the cone, the amount 
contributed is zero. The contributions from the plane 
waves starting out from the lifting surface will be sepa- 
rately considered. Integration over the wing chord there- 
fore gives for the downwash velocity induced by the lifting 
surface 

2 tt c z = 7 0 r f(i - S'.-n.Od v 
'-A ' = o 

or written out in full 
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The uppeT integration limit f 1 • = £ x . is different 
according to whether the point lies within the 
Mach cone II arising from the end point of the 
trailing edge (fig-. 5) or between the latter and 
cone I arising from the end point of the leading 
edge. The corresponding limits will he 

| x = 2t/b = € (within cone II) V 

l 1 = | - K 7(n - 1) S + f (between cones 1^ (31 ^ 

and II) I 

as may he easily seen after some consideration. 
Introducing the new variables of integration 

■ a - r) 2 - * t 2 = t 

and writing for briefness 

ai 3 = (T) - l) c + i* 

we have 



2 TT — — a — ^ ~ 



T = T 3 
- 1 /* ^ T - K& 

1 S 'J T W- K 



d T' 



T = T, 



T = T, 



(ti - I) 2 



]T_- 1 
a 



> 2 J 



a t 



K 2 (TI - 1 V 



T = T. 



T = T. 



d t 



(ti -1)' 



i 

O 

II 

'# 

03 
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The evaluation of the integral gives 



c z 

2 tt — = - 



7 o (r\ - l) a +- t 



' ■ • •• • ' v- h 



r. /e - "'I')' - «c B {<n- i)" + i 

- k| arc tan * ■ 



e»= o. 



Taking account of the different upper integration limits 
according to whether the point considered is within cone II 
or between cones I and II, equation (31), and setting for 
"briefness 

Ul € = (| - €) 2 - K 2 [(T1 - l) 2 + £ 2 ] (32) 

there is obtained as the final expression for c* z : 
For cone II 

2 tt !i = ILli J'vG - vG7) 

V <n - D a + t al 

- k j arc tan arc tan > (33a) 



Between cones I and II 

c_ r> - 1 /— VuT 

2 tt — = Viu - K arc tan (33b] 

y o (T) _ D 2 + t S «c<H - 1) 

where 

K ( T) — 1) _ . _ 

_ 1< _U — < + l . _ tt < arc tan < tt 

From these formulas it follows that <T on the surface of 

z 

cone I is equal to zero and on account of u; = 0 is con- 
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tinuous in passing through cone II. On the common axis of 
cones I and II, there still occurs the same singularity as 
in the case of the lifting line, cT z , there becoming in- 
finite as l/r • 

In order to obtain the total downwash velocity, there 
is still to be added to equations (33a) and (33b) the por- 
tion contributed by the plane wave. This contribution is 
different from zero only between" the plane waves starting 
from the leading and trailing edges (fig. 4) . According t 
equation (20) 



'wave 



2 y °' 



7wave 



= Ctt = 0 

x wav e 



(34) 



The expressions for the two remaining components of 
the induced velocity cf y and cT x are found by similar in 



t egrations. 
Cone II; 

2 TT — — =5 .— 



ffe shall only indicate the results: 

t 



On - D a + V 



2 tt — = arc tan 



<e - e)cn - i) 



- arc tan 



Between cones Z and II: 



t(r\ V 1) 



>(35) 



2 TT ~ = - 



2 tt = - arc tan 



(T) " 1) + V 

Mr) - 1) 




l/r) still remains, however, and prevents the computation 
of the induced drag for this' lift distribution. 
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Tor the downwash distribution in the plane z = O t 
at 1) the location of the wing x < t, and at 2) behind 
the wing x > t , taking account of the >plane wave, there 
is obtained the following: 



1) For x < t: 
2c 



'Z 



- 1< *<..(>:■ 



o "I 1 l & -s/i - $ 



= - 1 + 



- arc tan 



} 



2c. 



0< + 1: 



k7. 



— arc tan 



AT 



•} 



S (36a) , 



where for the arc tan the same values are to be taken as 
in (33) and. $ is given by equation (15a). 

2) For x > t: The plane waves do not contribute 
anything but the formulas obtained differ according as the 
region considered is within cone II or between cones I and 
II (fig. 5). 



. (i . *.)<«<! - i.^a- 

0 



■1 / Jl - •d 2 
- i ( arc tan arc tan 

tt \ 6 



) 



-><*<- (l-f) 

J 



- arc tan 



) 



(36b) 
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The downwash di str ibut<ion- fc f or x < t- and for x = 2 t 
computed "by the above equations is shown in figure 6. 

Further, We have in the same manner as for- the lift- 
ing line determined for the lifting surface the stream- 
line field of the y and z velocities in a plane at 
right angles to the cone axis* At the ' location of the 
wing (x < t, fig. 7) there is obtained outside of the 
cones springing from, the wing tips a constant downwash f 
due to the plane waves, along the span. The streamline 
picture within the Mach cone in the outer half is similar 
to that of the . lifting line (fig. 2); the inner half how- 
ever is entirely changed by the additional downwash ve- 
locity from the plane wave. 

The streamline picture behind the wing (x = 2t , fig. 
8) has, outside the Mach cones springing from the wing 
tips, a constant downwash velocity due' to the plane waves 
in two strips symmetrical- to the plane z = 0.' These two 
strips are limited by the plane waves starting out from 
the forward and the trailing edges of the lifting surface. 
Within the Mach cone the streamline picture in the outer 
ring is the same as for x< t and is changed only in the 
inner region. 

We shall yet consider briefly the question, what the 
form of the wing surface must be that corresponds to the 
assumed lift distribution. The wing plan form we have as- 
sumed as rectangular. Angle of attack and twist are ob- 
tained from the consideration that at the wing, i.e., in 
the plane z = 0 in each section parallel to the flow di- 
rection, the direction of the streamlines must be parallel 
to the wing tangent. Let z = z(x, y) be the equation of 
the wing surface and z(0, y) = 0, i.e., straight leading 
ei?e. Then we have 

d z c 2o (x, y) 
d x u 

o 

where c. includes the induced velocities from both the 
z o 

plane waves and the edge cones. There is thus obtained 
for the wing surface 

3C f = X 
1 /"* 

z(x f y) = — / c Zo (x«, y) d x» (37) 
u o J 
x' = 0 
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so that a further quadrature is required to compute the 
form of surface wing. 

For the case considered of constant lift distribution 
there is obtained for the region outside of the two Mach 
cones at the wing tips, from equations (37) and (20): . 

z(x, y) = -P 0 x 

that is, a flat surface with angle of attack p Q . Within 
the Mach cone the surface "bends downward more and more, 
strongly as the edge is approached. The edge itself (y = 
±b) is bent infinitely downward, i.e., actually the rec- 
tangular surface with constant spanwise and chordwise lift 
distribution is not possible. For this reason we may dis- 
pense with the further computation of the wing-surface 
shape. 

V f TRAPEZOIDAL WING WITH CONSTANT LIFT DISTRIBUTION 



We consider now a trapezoidal wing with constant sur- 
face density of the lift 7 Q (fig. 9). If the wing is 

cut away behind (taper angle t. fig. 9) in such a manner 
that the Mach cone at the tip of the leading edge does not 
overlap the wing (t > a) , the induced drag is obviously 
equal to zero and only the wave drag exists (reference 4). 

t> a:Wi = 0 

The trapezoidal wing with constant surface density of the 
lift 7 0 is plane outside the Mach cone and has the angle 

of attack p o where 

7 o = 2 f*o u o tan a 

The trapezoidal flat surface with constant lift distribu- 
tion whose cut-away angle t is greater than the Mach 
" cone angle may be looked upon as the~ "ideal supersonic 
wing with finite span" sin,ce for it the ratio of drag to 
lift is no greater than for the w}ng of infinite span. 

The computation of the induced drag for t < a is 
possible in a simple manner from the above results. By a 
lifting element we shall mean a strip of the lifting sur- 
face of chord d x and therefore with circulation 7 0 d x f 
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Such a lifting element at x » 0. generates at a lifting 
element of chord d x 1 at x = x 1 

, y=yi(x) 

42 ff iox« = p 7 o 4 * f 4 c i° 0 3C,) 4 y < 38 > 

y=y 0 00 

(ox') 

where d c ZQ denotes the downwash velocity induced hy 

the lifting element x =.0 at the position x - x 1 . The 
integration limits are the surface of the Mach cone aris- 
ing from the tip. of the wing leading edge and the side 
edge of the plate. ?or the downwash velocity 
e Cox») 
zo 

(14c) 



in the plane £ = 0, we have according to equation 



(ox-) = 2sl 7e 3 - ifOn 



d 'XT' = "T 1 S " * * (39) 

with the aid of ..which equation (38) , becomes 



d= i = -£-V ««*■«' /' -Z^-" 1,2(11 - x) - d „ 



2 n " ?(n - i) 



'0 



T| - 1 

or with d =K , according to equation (15a) and 



i 



e = (40) 

tan a 



as the reduced angle of taper 



d *iox' = - — 7 Q 3 d x d x' / - ' T - dd 

2 TT • J d 

d=6 



d=i _____ 



-£- 7 0 3 d x d x' g (6) (41) 

c TT 
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The evaluation of the definite integral gives 

g (9) = - fl - 6* - log 1- A - 6« (42) 

■ ■ e 

According to equation (41) the induced drag from the lift- 
ing element x = 0 at the position x = x ! is independ- 
ent of the distance "between the two elements. All ele- 
ments lying between x = 0 and x = x 1 accordingly pro- 
duce the same drag, so that the total drag induced at x = 
x 1 amounts to 

d W ix » = - JL 7 Q 2 x' d x' g (6) 

2 TT 

The drag for the entire wing is obtained from the above 
"by integrating over x 1 "between the limits x 1 = 0 and 
x« = t and multiplying "by two (both ends) 

x f =t 

Wi = - £. 7 0 2 g (6) j' x' d x» 

x'=0 

= - JL 7 0 2 t 2 g (6) = - JL r 8 g (6) (43) 
2 tt ° 2 tt 0 

The minus sign is explained by the faict'that with our 
choice of coordinate system the drag component of a force 
is in the direction of the negative z axis. Formula (43) 
for the induced drag of a trapezoidal wing with constant, 
surface density of the lift is of the same structural form 
that is found for the incompressible flow. For triangular 
lift distribution (lifting line) in the case of incompres- 
sible f low, we have , _f or„exampl_e , _ 

log 2 „ 2 

TT 

where T 0 is the circulation at the wing center. For 

equal total circulation r o , ffj according to equation 

(43) is independent of 6, i.e., the ratio of the tangent 
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of the angle - t to the tangent of the Mach angle (equa- 
tion (40)). In passing to the rectangular wing, 8 > 0 f 
the induced drag according to equations (42) and (43) he- 
comes logarithmically infinite; in agreement with our re- 
sults of the previous section. 



Actually, we are not interested so much in the value 

of the induced drag alone as in the sum of the induced and 

wave drags. For the wave drag, according to equation (21)* 
we have ' ' 

*wave " p F Y o c z* avQ 



where F = h t 



hence 



t tan t \ 



the area of the wing 



> (44) 



'wave 



wave 



2 tan a 



P 
2 



y 0 3 

tan a 



For the lift we have, on account of 7 Q = 2 p Q u Q tan a: 



A=P7 0 u 0 F=2P u 0 2 F p 0 tan a 



(45) 



or 



£ u 0 * 

2 0 



= c a = 4 p 0 tan a 



(46) 



For the wave drag we obtain from, (44) 

• "wave = 2 PI u 0 2 p^ 2 tan a 



'wave 



= c. 



'wave 



= 4 p 0 a tan a = p 0 c a 



(47) 



and for the Induced drag from equation (43) 
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ffi = JL t 2 4 p Q 2 u Q 2 tan 8 - a g (6) 

2 TT 



*i _ 1 4 t j3 0 2 tan 2 a g (6) 

— — — c, 



£ u 8 y rr ■ /. . t tan t 



71 D (i . t tan T ^ 



c W4 = 4 Po 2 tan a - JL<|L (48) 
w i 0 tt i _ 6 X 



where X = * ■ * &n °^ is the "reduced aspect ratio" of the 
b 

wing. For the total drag there is thus obtained from (47) 
and (48) 

<cw) wa ve + ind = 4 Po" *« * {* + " r^^TT } 

TT 1 - O A 

/ x c a 8 - X g (8) ^ 

< c w) W ave + ind = U + } < 49 > 



It follows therefore from the above that for supersonic 
speed the wave plus induced drag, like the induced drag 
in the incompressible flow case is proportional to the 
square of the lift. Equation (49) is analogous to the 

C 2 p 

well-known formula c w . = — ^ of the elliptic lift ■ 

1 tt b 

distribution for the incompressible flow. The essential 
difference lies in the fact that for the supersonic flow 
the drag parabola for small aspect ratios t/b is to a 
first approximation independent of the aspect ratio. The 
manner in which the drag increases" wi th increasing reduced 

aspect ratio X = » ■ tan °^ and decreasing 8 is shown in 

b 

/ 

figure 10 where c w / is plotted against A for 

/ 4 tan a 

various values of 8. Our formulas are valid only for 
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X< — , i.e., for the case in which the Mach cones do not 
2 

overlap on the wing. 

In order to be able to predict what the wing shape 
must he sc that our assumed lift diBtribut i on may he pos- 
sible, we must first compute, the field of the induced 
velocities. For this purpose equation (39) is to he in- 
tegrated over the trapezoidal area. The value 

ic l^^ according to (39) gives the downwash velocity in- 
z o 

duced at the position £ by a lifting element 7 0 dx 
starting at £ = 0 and ending at T| = 1. A lifting ele- 
ment which starts at £ = £ 1 and ends at r\ = r\* thus 
produces at the position r\ f £ = 0 the downwash veloc- 
ity 



z o ^ 

.2 tt (5 - *•)(*.•■ - r\) 

For the velocity induced by the entire surface there is 
thus obtained 



4*> = - 2sl / ^TT^T') 2 - « a 5 - ^i') 2 4e . (50) 

0 2 71 J : {% - | 1 ) (*n* - tj) 

I '= o 



In order to evaluate this integral we introduce the new 
integration variables 

* =K T D | 3 K tans)'' (51) 

£ - * 

(8ee fig. 9.) Since the end points of the lifting elements 
lie on the wing contour there exists the relation 

' . - • • 

1 . n» = £' tan T (52) 



The upper integration limit £ ' = £j, in equation (50) is 
obtained from the condition. (See fig. 9.) 
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tan = tan a: | ! = : d f = 1 

The lifting elements whose ? ! is greater than the l x 
thus determined give no contribution at the points £, T|, 
I = 0. For the lower integration limit. 

£« = 0 : *n« = 1 : j^K I! =r * 

From equations (51) and (52) there is obtained 

d g I d fl 1 

e - e 1 " « f - 6 

where 8 is the abbreviation introduced in equation (40). 
There is then obtained from equation (50) 

, ( d ) = - 2l k /' ^ x I *'* d = k j( di e) (53) 

In evaluating the above integral the following three cases 
are to be distinguished: 

1. 0 < 6 < a ; 2. 0 < $ < 6 ; 3- d < 0 < 6 

In case 1 the point P(fl) lies within, in cases 2 and 3 
without the trapezoidal wing. . In case 1 the integrand is 
regular over the entire range of integration; in case 2 
it possesses a singularity at d 1 = &; and in case 3, two 
singularities at « 1 = 0 and $ ! = £. In cases 2 and 3 
the principal values are to be taken, namely, 

0 < * < 9: tf'^e-c j 6 f = i 
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a < 0 < 6: 
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fU 6) slim 4 / d d' + 

- v J *•(#• - 6) 



e— »t>; 



The integral 



♦X* 1 - 6) 



?.(*.e) = y „ . d (55) 



may be obtained by elementary methods. We set 



2 = t tf« - 1 



where t is the new integration variable so that (55) be- 
comes 

t = t 

... (1 - t 3 )2 dt. 

FU.6) = 



t =i 

where *. 



J t(l + t 2 ) {2 t - 6 (1 + t 3 ) ) 



1 + J 1 - a 2 /.-n 

tj = v|/ = — (56) 



By breaking up into partial fractions there is obtained 
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(1 - t 8 ) 



2\ 3 



t(l + t 2 ) [2 t - 6(1 + t 2 ) j 



1 + t 3 tt 



v t - t B t - t/ 



where 



1 ± 



yr: 



■ 2 



1,0 6 
Performing the integration, there is obtained 



- 2 arc tan t - 



log 



t . J\ - 6 2 



1 t = 

{log (t - t 2 ) - log ( t - tj) J J 



1 



For 0 < 6 < $ there is therefore obtained directly 



( dt6 ) =E-2 arc tan *- 

2 6 



777 



log 



^1 + Ji - e 2 v+/ e - 1 + /1 6 C 



(57) 



1-^6 + /i - e 2 



while the formation of the principal value according to 
equation (54a, h) gives 
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0 < d < 6: 



log ^ 



\ 



*U,6) = ! - 2 arc tan + - ±2| 



7l - 6 a A +Jl - 6 a <|r 9 - 1 +/l - & 2 
log 



+ 8-l-</l- 6 s 



* < 0 < 6 : 



\(57) 



,6) = - E - 2 arc tan + - 
2 



log (-+) 



_ 6 8 • /i - 6 a ^6-i +71 - e a 

+ log . 

6 \ 6 *e-i-/i-e«/^ 



where 



- H < arc tan V < II 



There is thus found the downwash distribution in the entire 
Mach cone springing from y = b/2, x = 0. For F«,b) we 
have 

d = ±X:?(*L; 6)5 0 

d = 8, ^ = 0:F(8,6) = F(0. 8) = » as log a at d = 0 (58) 

On the two rims of the cone (d = *D the induced velocity 
is thus zero and on the edge of the trapezoidal wing 

'= 8) and on the cone axis U = 0) it is infinite. 
In figures 11 and 12 for the particular case 6 = 

i (tan a = Jz, tan j = — ) there is shown the induced 

3 v J~z' 

downwash velocity in a section parallel and perpendicular, 
respectively, to the principal stream direction. 

To the above velocity field of the tip vortices there 
is still to be added the velocity field due to the plane 
wave. The latter in the plane z = 0 within the wing area 
is 



N.A.C.A.. Technical Memorandum No. 897 29 



= - i -^ L - = -eo*o 

° wave 2 tan a 
and outside the wing area 

c z°wave B 0 

Prom the velocity field it is now possible to compute 
the form of the trapezoidal wing surface that has constant 
lift distribution. Outside of the Mach cone we have, ac- 
cording to equation (37) 

x 1 = x 
x' = 0 

that is, a flat surface with the angle of attack {3 Q giv- 
en by equation (18). The twist of this flat surface within 
the edge region of the trapezoidal area that is overlapped 
by the Mach cone is given by 

z(x, y) = ~ J c 20i d x' 
x'=(b/2-y) /tan a 

and according to equation (53). 

x'=x x'=x 

z(x,y) = -i- 2l K I F(*« ,6) d .x> = -?i /* F(0«. 6 ) d X' 

2TT U Q J TT J 

x' = (b/2-y)/tan a x' = (b/2-y) /tan a 



On account of 3 = - there is obtained 



'=0 

?($■ , e) 



^18 



d V (59) 
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Since the function 6) is known from equation (57), 

it is possible from the equation above to compute for a 
given 9 the profile sections of the surface at various 
distances (b/2 - y) from the edge. The ordinate of the 
obliquely cut-away edge of the trapezoidal 'area for b/2- 
y < t tan t: 



'=9 



( , B ,, )= ^»(|. y ) ; y-^ 44 . 



(60) 

d' = i 



The integrand becomes infinite for = 6 (equation (58)). 

The integral exists, however, and may be evaluated by spe- 
cial, computation. There is obtained 

a'=6 



(60a) 



/' 9) j A t _ _JT arc sin 6 l_ .1 

d'=i 

(The evaluation of the integral was performed by Dr. P. 
Riegels.) 

For 6 = l/3, we thus have 

6 = 1/3: =-4.92 

The ordinate of the rear edge point x = t , & - y = f tan t 

for 6 =■ 1/3 is thus 2 = - 1.532. £ d t. (?lat surface 
z = -p Q t, twist z = -0.522 p 0 t.) 

For the special case G = 1/3 (tan a = *fz ; tan T = 
— ~ ) the profile sections have been computed and are 

given in figure 13. If the trapezoidal wing were flat 
there would be a drop of the lift toward the edge down to 
zero. In order that full lift be maintained up to-, the 
edge, the wing must be bent downward. The twist of the 
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wing 
from 



directly at the edge is very strong as may he seen 
the "elevation contour lines" (fig. 14). 



VI. COMPUTATION OF THE LIFT DISTRIBUTION 



FOR THE UNTWISTED RECTANGULAR WING 



The examples thus far considered are all in connec- 
tion with the so-called first principal problenT of the 
airfoil theory where the lift distribution is given and 
it is required to find the drag and the wing shape. Of 
greater practical importance is the second principal prob- 
lem where the wing shape being given it is required to 
find the lift distribution and the drag. As in the case 
of the incompressible flow, eo also in the case of the 
compressible flow the first problem, which leads only to 
quadratures, is considerably more simple than the second, 
which requires the solution of an integral equation. 

In what follows there will now be given an example 
of the second principal problem, namely, the computation 
of the lift distribution for a plane rectangular wing 
(span = b, chord = t), that is to say, the same problem 
that was first considered by A. Betz (reference 5) for the 
case of incompressible flow. In the treatment of this 
problem we can utilize to a large extent the results we 
had obtained in the previous section for the trapezoidal 
wing with constant surface density of the lift. We con- 
sider a rectangular flat plate which extends from x = 0 
to x = t and from y = -b/2 to y = +b/2 and is set at 
the small angle of attack p 0 to the undisturbed veloc- 
ity u Q (fig. 5). Within the region bounded by the. plane 
waves starting out from the leading and trailing edges and 
the two Mach cones there is the constant downwash velocity 
due to the plane waves 



Outside the region of the flat plate overlapped by the. 
Mach cones at the tips there thus exists the constant lift 
distribution 7 0 . At the tips y = *b/2 the lift must 
vanish, that is, 7=0 at y = ±b/2. There is required 
the lift distribution 7=7 (x^y) within the region 

overlapped by the Mach cones. The problem is considerably 



_ c 



zo wave * 




(.61). _ 



32 



N.A.C.A. Technical Memorandum No. 397 



simplified by the circumstance that, as will immediately , 
become apparent, 7 does not depend on. the t^o independ- 
ent variables x, y t but only on one of the variables 

b 

*--2 (51) 

. x tan a 

(fig. 11), For the required lift distribution 

7U) = 7 Q f (*) (62) 

of the rectangular wing there then exist the boundary con- 
ditions 

, 4=0;f(i) = 0 1 

> (53) 
0=1: f («) = 1 J 

In order to be: able to set up. the integral equation for 
7(^) we must first compute the field of the downwash ve- 
locities w(#) induced by a. rectangular wing with the cir- 
culation distribution 7(£) in. the plane z = 0. The in- 
tegral equation for 7(#) is then obtained, in the known 
manner from the consideration that for each position of 
the wing the sum of the effective angle of attack 

- i tC*) , : 

^~ (54) 

2 u 0 tan a 

and the induced angle of attack - w Lil must be equal 
to the geometrical angle of attack £ 0 

- = p 0 (65) 

The velocity field w(^) induced by the edge vortices is 
obtained by considering the rectangular wing with the var- 
iable lift distribution 7(d) = 7 0 f U) as built up by the 
superposition of trapezoidal wings with various taper an- 
gles each of which wings possesses a constant lift distri- 
bution.^ Again, let 6 = tan < T be the "reduced taper angle 11 

tan a 

(equation 40), then the lift distribution 7 = 7 Q f (d) may 
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"be obtained by the superposition of trapezoids with angles 
6 and lift densities 7 0 f'(6) d 6- Each of these trape- 
zoids produces, according to equation (53) the velocity 
field 

a w ( d > = „ 7 ° fl (6) I (i,6) d 6 
2 tt tan a 

and integration over 6 from 6 = 0 to 6 = 1 then 
gives the induced velocity field over the rectangular wing 

6=i 

w U) = 7 ° /" f'(6) P U,6) d 6 (66) 

a J 



2 tt tan 

6= p 

By substituting the above expression for w(£) in equation 
(65), there is finally obtained, taking account of (61) and 
(64) the required integral equation for f(#)s 

6=i 

f(«) + J y f 1 (6) F («,6) d 6 = 1 (67) 
# 6=o 



to which are added the boundary conditions (53). This 
integral equation for the lift distribution has the same 
structural formas that for the incompressible flow. It 
differs from the latter, however, by the different core 
F($,6), which is given by equation (57), and for the super- 
sonic flow is of a much more complicated form that for the 
incompressible flow. Equation (67) also exhibits the nota- 
ble property that neither the aspect ratio of the wing nor 
the Mach number appears explicitly, whereas in the incom- 
pressible case the characteristic value of the integral 
equation depends on the aspect ratio. ThedepencLenceof 
the lift distribution on the Mach number appears in the 
introduction instead of the geometric angle <p (fig. 9) 

tan qp 

the reduced angle # = as the variable. It is neces- 

tan a 

sary to solve the integral equation (67) only once to ob- 
tain the lift distribution of the rectangular wing for all 
aspect ratios and all Mach numbers. 
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The solution- of the integral equation (57) appears 
at first sight .quite diff icult ,. particularly . on account 
of the complicated structure of the core 6). (See 

equations 56 and 57.) By a simple transformation of 
equation (67) it is possible, however, to simplify the 
problem considerably.* The equation is [a. nonhomogeneous 
integroiif f erential equation for f(#). Instead of it 
we shall consider the equivalent equation for f'(d). 
Taking account of the singularity of the core, equation 
(67) may be written 

6«* 6=i 
fU) + -{ J f'(e) FU, 8) d 6 + j f'U) 1(4,8) d 8 } = 

6=o 6 = * 

Differentiation with respect to d gives 

• 6=3 

f •(*)+ i ■{ f ■(*)*(*.*) + / f'(6)iid 6 

TT , J d * 

8=o • 
6=i 

- f'U) 3) + /' f(6) — d e}= 0 

J d .3 J 

and because 

d P. • Jl - 3 s 



d 3 «( * - 6) 

according to equation (53): 

— - . ' 

„(.) ..i^ii*! /' d 8 = 0 (68) 

e= 0 



♦For this suggestion I" am indebted to Doctor Lotz and for 
carrying out the numerical solution of the integral equa- 
tion to Mr. Pret8ch. 
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The above is the equivalent int egral equation for f'(d) 
which, however, is now homogeneous . : The solution of this 
integral equation for f'(6) is possible by building up 
f 1 ( £ ) in n steps and solving the corresponding system 
of linear equations 

f '<W - . ' • ' 



2V+1 



17 




= 0 

(V=0, 1, n - 1) 

(69) 

This is a system of n homogeneous equations for the n 
unknowns f 1 ( £ 2 u+j) ( V = 0, 1 , . . . , n - 1) . Since, as 
closer investigation shows, f'(0) = f f (* ) is suit- 
ably chosen not constant but equal to 

f t = a - b ^ 

There is then obtained in place of equation (69) a 
nonhomogeneous system of equations of the nth order for 

the n unknowns ^ , (i 2 y +l ) (v = 1 . . . , n - 1) . The 

a a 

further unknown a is obtained in the numerical integra- 
tion for f ($) from the condition 



u=o a 



(70) 



In carrying out the numerical process there were first 
taken five steps (tf 2 v+i = °.l; 0.3; 0.5; 0.7; 0.9), then 
ten steps C* a u+x = °- 05 5 "0.15; . ..; 0.95). It was found 
that the ten-step approximation gives an improvement over 
the five-step process only in the interval 0 < 3 < 0.2. 
In the third approximation therefore only the interval 
0 <.*<0.2 was again subdivided (^ 3V+1 - 0.025; 0.075; 
0.125; 0.175). The values obtained in thi s manner for 
f'(£) and f(fl) are given in table 21 and the function 
f U) plotted in figure 15. At $ = 0 the function 
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f (#) * possesses a .singularity since ^ f. 1 (d) there becomes 
infinite. The mathematical nature: of this singularity 
could not as yet he determined* 

V7e shall now compute the lift, wave drag , and ..induced 
drag as well as • the moment about the transverse axis of 
the rectangular flat surface. 

The lift Aj of that portion of the surface which 
lies outside the two Mach cones is 



/ t tan a \ 
Aj = P u Q 7 a > t (l ) 



b 

while the lift of the two triangular portions overlapped 
by the Mach cones is 



A II = p u o t3 tan oc J 
where 



'V A A — n „ 

0 'O 



7 d d.= p u D 7 Q t a tan a Z 



£ = / f («) d i = 0.684 ,. (71) 



i) = 0 

The total lift of the rectangular plate is therefore 

A = P u 0 7 Q D t { l - (1 - K) x} 

or, according to equation (20) 

A=2pu 0 2 j 0 nan a {l - (1 - X) x} (72) 

For the lift coefficient there is thus obtained 

c a = 4 p 0 tan a {l - (1 (73) 

For ;the -wave drag'outside of' the Mach cones there is 
obtained simply • -V 



p 2 - b t -t tan a ^ . v 
wave I = ?o A i = | V {'I } < 74 > 



tan a 
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The wave drag of the two triangular portions overlapped 
fcy the Mach cones is 

w wave II = 2 P 1 7 c wave d f 



where 



' 2 wave 



f («) 



2 tan a 2 tan a 



and 



d f = - t s d (tan cp) 
2 



Table II 

Lift Distribution of the Untwisted Rectangular Wing 

fU) and f '(d) 





f •(«) 


6 


f (*) 


0 


00 


0 


0 


0.025 


4 .49 


0.05 


0.219 


.075 


1.86 


.1 


.312 


.125 


1.39 


.15 


.381 


.175 


1.24 


.2 


.444 


.25 


1.10 


.3 


.. 554 


.35 


.958 


.4 


.649 


.45 


. . 850 


. .5 . 


.734 


.55 - 


.753 


.6 


.810 


. 65 


. 655 


.7 


.875 


.75 


,546 


.8 


.930 


.85 


.417 


.9 


.971 


.95 


.225 


1.0 


1 



We then have 



'wave II 



= - — !L- /W 8 d (tan cp) 
2 tan a J 



t* 7 0 3 

2 



P 
2 



'wave II = * V *' £ i 



f 2 d a 



) = 0 



(75) 
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wnere < - - <S - i 



d = o 

Similarly there is obtained for the induced drag in 
the two triangular .regions overlapped by the Mach cones 

^ = 2 P / 7 c 2l d f 

where from equations (62) and (65) 

c zi = i-^- (i - fU)> 
A 2 tan a 

We then have 

" P " H ' 

3 v 2 > fU) [l - fU)J d« 



r t = - t a 7 0 a (K - K x ) (76) 



For the total drag 

W = w wave I + ? wave II + w i 
there is thus obtained from equations (74), (75), and (76). 



P a b t 
I - - 7 0 a 



2 tan 

or from equation (20) 

= 2 P u 0 3 |3 0 a P tan .a |l - (1 - K) X } (77) 

and for the drag coefficient 

c w = 4 p 0 s tan a {l - (1 - K) X j (78) 

Prom equations (72) and (77) there is obtained between the 
lift and the drag the simple relation 

W = j3 0 A (79) 
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There is thus obtained for the plane surface of 
finite span the same simple result as for the infinitely 
long flat plate, namely, that the ratio of the total drag 
for a frictionless flow to the lift is p 0 : 1. This may 
also he explained by the fact that in contrast to the in- 
compressible flow no suction force arises at the leading 
edge in the supersonic case and the resultant air force 
is therefore at right angles to the plate. 

For the relation between the drag and lift coeffi- 
cients, there is obtained finally from equations (73) and 
(78) 



c w = 



'a 2 1 c a 2 



4 tan a 1 - (1 - K) X 4 tan a 1 - 0.316 X 



(80) 



The above formula has the same structural form as 
formula (49) for the trapezoidal wing with constant lift 



distribution. In figure 10 c w / = has been plotted 

/ 4 tan a 

against the reduced aspect ratio X (dotted curve). It 
may be seen that the rectangular plane wing for the same 
lift has the same drag as the trapezoidal wing with con- 
stant lift distribution with the reduced taper angle 6 = 
t an t 

= 0.27. . For the reduced aspect ratio X= 0.3 the 

tan a 

rectangular plane wing has, for the same lift, about 10 
percent and for X= 0.5, 19 percent more drag than the 
ideal trapezoidal wing whose taper angle is greater than 
the Mach angle. 

With the above results the theoretical polar and 
moment curves for the plane rectangular wing may be given 
for various aspect ratios and Mach numbers. For the mo- 
ment Mg about the transverse axis in the wing leading 

edge , there is obtained- - _ . . . 

Mg = 2 P u 0 2 p 0 tan a b 1 2 { i - | (1 - K) X } 

* It is interesting to note that the constant 1 - K = 1 - 

J X t ( 6) d 3 is equal to 1/tt within the computational accu- 

0 /. 

racy. That this is exactly so has as yet not been shown. 
For this it would be necessary to know the exact solution 
of the integral equation (67). 
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and for the moment coefficient c 



»H 



'm 



H 



m. 



H 



= 4 Po tan a{|- £ (1 - K) X} 



. § u 0 3 d t* 



m H a 



I - 0.211 X 
1 - 0.316 X 



(31) 



Through equations (73), (80), and (81) f- the polar and mo- 
ment curves not considering the frictional drag, are com- 
pletely determined. In figure 16, the polars are given 



for the aspect ratios 
u 



1=0, i, and 4 and for the Mach 
b 5 2 



numbers -2.= 1.2, 1*5, 2.0, and 3.0. 
c 



The drag differences 



"between wings with various aspect ratios are considerably 
smaller in the case of the - super sonic flow than for the 
incompressible flow since in the first case the greatest 
part, of the drag is contributed by the wave resistance, 
which is independent of. the aspect ratio. 

The plane rectangular wing at supersoni c f low is one 
with constant center of pressure position, if the fric- 
tional drag is disregarded. The position of the center 
of pressure depends only to a slight extent on the re- 

t ian a 

duced aspect ratio X = 1 Tor the infinitely long 

b . 

wing, the center of pressure lies at the midchord position 
and with decreasing aspect ratio it moves forward somewhat 
(table III). 



Table III 



t tan a 

"b 


0 


. . . V 5 




° m H.. 


1 










0.489 


6.469 




I 
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Formula (80) for the rectangular flat plate is the 
analogy to the familiar . c^. =- c y 2 . F/tt b 2 of the incom- 
pressible flow. Like the latter it enables the recompu- 
ta^ion of the drag from one reduced aspect ratio 



t tan a, 



to another X 2 = 



tg tan as 



From equa- 



tions (73) and (80) there is obtained for the new angle 
of attack and the drag 



- -P. + t i 



tan a a (1-0.316 \ 2 ) tan' 04 (1-0.316 X t ) 



} 



c a 2 j 

;w 2 = 0 ^ + T V 



_J } 

tan a (1-0. 316 X 2 ) tan a x (1-0.316 *J 



!> (81a) 



VII. TRAPEZOIDAL LIFT DISTRIBUTION 

a) Lifting Line 

As a further example we now compute the induced drag 
and the velocity field for trapezoidal lift distribution, 
for both the lifting line and the lifting surface (fig. 
17). Let the lift distribution therefore he given hy 



ron') 1 - «n« 



1 - n : 



for Tlj £ «n« < 1 

_ 1 < T)« < - T| 1 



r(-n') = r for - , n 1 < r\ < + r\ l 



(82) 



where Tj. = b'/b, according to figure 17. The field of 
the -induced velocities and- induced drag for variable „lift_ 
distribution may be obtained in the familiar manner from 
the lift distribution by superposition. On account of 
integration difficulties, however, this computation can 
not directly be made on the potential but must be carried 
out separately for the three velocity components. From 
equation (14c) we have for the induced downwash velocity 
of a lifting line ending at T\ = TV with circulation r Q : 
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c z (T!») = -2. : , ,;, — (83) 

TTb (ja-Kaj»)[flrKn>X*|.J a l/ea.. lc a Rn'-'n) a + S a 3 

From the above there i.a obtained by superposition 
the downwash velocity "~ c z for variable circulation r(ty ) : 

Tl'rr+l 

c z .= - ~ /" o B W) L£ d Tl« (84) 

< . • 0 v=-x " 7*7: " 

'For the trapezoidal lift distribution according to 
equation (82) we have therefore if, on account of symmetry, 
.we restrict ourselves to the* half-wing y > .0 

c z = - ■ / • • c«(ti») d r\* 



1 - -V 



or, according to equation (83) 



b c z 

TT 



or with 

ks (n» - t;) 3 + t a = t 



* Cz /». d t 1/2 {; /> d T 

r o 1 - nrJ We 2 V 1 - * 8 - "£ 2 ^ V| 2 - T 

Performing the integration there is obtained for 

points £, t], I within the Mach cone at the wing. tip r\ = 1 



b c 2 _ 



r o 



.'-JL.{.-jjs5L + Ii. g i±i£'} (ss) 
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and a corresponding expression with reversed sign for the 
Mach cone- at T| = T} r ; The value of cu 1c. here given by 
equation (12), In the cone r\ = 1, c 2 > 0 so that there 
is upwash velocity. In the cone r\ = r\ 1 there is a 
downwash velocity of the same absolute magnitude (c z < 0) 
and outside of the. two cones c z s 0, a result which, is 
also to be expected from reasons of symmetry since, on ac- 

d r 

count of = cons t • t * all separating vortices are of 

d 

the same strengths With 



*(T1 - 1) 



and = 



K(T1 - TU) 



there is obtained for the downwash distribution in cone I 
and III respectively in the plane z = 0. 



TT 



*> Cz . 



i r r a 1 1 + vl-" A a \ 

= ± 1 \ - J\ - * + - log ■ . : \ 

' i - ^x 1 2 i - JT7T* 5 



(86) 



(86) 



On the cone surface according to equations . (85) and. 
c z s 0 jand is therefore continuous, in passing 
through the cone. ( On the cone axis c z now becomes log-' 
arithmically inf init e , . whereas with the rectangular .lift 
distribution (horseshoe vortex) c 2 becomes . inf init e pn 
the axis as. r" 1 . The • logari thmi c singularity, of c z is 

no longer a disturbing factor for the computation of the 
induced drag. 

.For the. sake of completeness there will also be given 
the remaining two components of the induced velocity. 
There is found for the cone at T| = 1: 



b c 



TT 



TT 




i - tu 



£(n - D 



(87) 



C 
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and corresponding expressions with reversed signs for the 
cone at T| s T| . For the arc tan there is to he taken the 
principal value 0 < arc tan < tt. For the outer cone 
(r\ as 1) the arc tan. is zero in ; the upper half plane on 
the outer quadrants, of. the cone surface and equal to +tt 
on the inner quadrants* In the wedge-shaped space be- 
tween the two cones c y is constant, being equal to 

p 

t <Of Cyrl— 2 ; Cx 2 0 (88) 

9 b - b' 

In passing through the plane £ = 0 t therefore 
there is a discontinuous increment in c v by 

2 . The region of the £ plane limited by the cone 

» - /■ b - b'x 

axes r\ = T| 1 and T) = 1 ^distance ■ j is thus a vor- . 

tex surface with constant circulation density the total circula- 
tion of which is equal to the circulation T 0 of the 
bound vortei in the region of the constant lift. 

A streamline picture of the y and z velocity 
components for a plane x = constant that intersects both 
cones is drawn in figure 18. Like the streamline picture 
for the constant lift distribution (fig. 2) it was ob- 
tained by computing the field of isoclines. On the outer 
halves of the cone surfaces c y and c z are equal to 
zero but the directions of the streamlines c z /cy have a 

value different from zero. In this case, too, not all 
streamlines are closed, part of the streamlines entering 
from the undisturbed region into the one cone and coming 
out from the other again into the undisturbed region. 



b) Lifting Surface 

In order to compute the induced drag for the trape- 
zoid-shaped lif t . di._stributi.on, we must,, as i» jwctlon IV, 
make the transition from the lifting line to the lifting 
surface. A rectangular lifting surface will therefore 
now be assumed of span b and extending from x = 0 to 
x = t. The chordwise circulation distribution is assumed 
to be constant of density T/t, while along the span the 
distribution is that given by equation (82). For the com 
putation we may here restrict our selves to the region be- 
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tween the cones springing from the leading and trailing 
edges of the lifting surface, since only this region en- 
ters into the question of the computation of the induced 
drag. We likewise need carry out the computation only 
for the cone at T\ = 1 ; for the downwash in cone r\ = T) 
there is obtained the corresponding expression with re- 
versed sigu. 

For the induced z component c z of the lifting sur- 
face, there is found, according to equation (85), with 
= r oA: 



e'-tx 



c z i ■ r l ( t „ V )J(^o^K[(^ ])«+ in _ 

c TT = • / d I 



1/2 /•* LzJll-ALz v) s -«* [ 3 t , 

+ / log r— d t 

i - tu j f. - e' - </( e - e«) a - * 2 C...3 
e'=o 



where 

according to equation (31). With the new integration var- 
iables I - £ f = I*, the above equation becomes 



7 o 1 - \ J e* 2 - * 2 £ 2 

. . _ 1/2 ..;£*. . . 

/ ^io g u* +y?* 2 - « s c — 3) d e* 

1/2 ' />£i* ; : 

log (t* - Ve* 2 - k 3 [- • • • 3) d i* 

where 



+ 1/2 /° • 
1 - 4*=i 
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• | x * = K »/ (T| - l) 3 + I* 

The three integrals are evaluated as follows: 
Setting f» a - K a £ a = T, we have for J x 

t 3 / ■ 

j = */ 2 /' v T — t 3 ( , n - i) a 4 T 

j; = ' i-/uT- K (Ti - 1) arc tan V 



With 



and 



— ■ lVUI • \ ' I — / — - ^» " — — 

1 _ r\ K . kC'H - 1) 



(t, - l) a + £ a = a^ 



there is obtained for J a 



1 1 7» a T a - K a a, a 
j = / i- log T d T 



x 

After a brief intermediate computation we have 
1/2 



|kE(1 - t,) 8 + i s 2 log (k7(1 - T!) 2 + C«) 

1 - *n i • 

... - I log (| +VuT) + Vuj. I 



and similarly 

1/2 



J 3 = + 



{ft[(i - -n) 3 + £ a ] log (k V(i - t}) s + c 5 -) 

- i log (e - JZ) + /^T| 
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By adding we obtain 

2 tt — = J- Vo) + k(T) - 1) arc tan 21 

1 - v\ x L ; K(Tl - 1) 

+ - I log — ) (89) 

{ - VU) 

A corresponding expression with opposite sign is obtained 
for the cone T) = T| . The arc tan in equation (89) lies 

within the range - — < arc tan < + — as follows from the 

2 2 
fact that c z must he symmetrical in (r\ - l) since the 

same holds for c z according to equation (85). 

The induced, z component thus found for the rectangu- 
lar lifting surface with trapezoidal lift distribution 
has the same singularities as the corresponding formula 
(85) for the lifting line. On the cone surface c 2 s 0 

and on the cone axis logarithmically infinite. For the 
downwash distribution at the location of the wing in the 
plane £ = 0, there is obtained 



o t j / 1 -/l - * 2 

2 tt =5 ± -VI - + $ arc tan + 



+ - log — — — \ = g ( d ) (go) 

K(T| - 1) K (Tl — Tl, ) 
where # = for cone I and # = i— for 



cone III, the upper sign holding for: cone I and the lower 
for cone III. Equations (89) and (90) include only the 
downwash velocity induced by the edge vortices. In order 
to obtain the field. of the total downwash motion, there is 
still to be added induced downwash velocity due to the 

plane wave. In the wedge-shaped space between the leading 
and forward edges of the wing (fig. 4) f this induced veloc- 
ity component is 
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'zo 



wave 



1/2 



2 t 



1 - n 



i - •n 1 



For the total downwash velocity in the plane z = 0, there 
is thus obtained from equations (90) and (91) 



For cone. I j 



- 1 < a < 0: 

0 < d < + 1: 

For cone III: 

- 1 < a < : 0: 

0 < d < + l: 



2(1 - tu) _ 



TT 



= e 



TT 



> 



(92) 



2(1 - -Hi) -. 
■ c 



zo 



^ TT 



The downwash distribution thus computed is plotted in 
figure 19. 



We are now in a position to compute, 
with trapezoidal lift distribution, the i 
order to avoid special complications, we 
the Mach cone springing from. the leading, 
does not extend beyond the wing tip and d 
the region of dropping circulation of the 
The first is identical with the condition 
springing from T\ = 1 does not extend in 
the wing whero the circulation is constan 
for the Mach angle the two conditions 



for the wing 
nduced drag. In 
shall assume that 
edge at r\ = r\ x 
oes not overlap 

other half-wing. 

that the cone 
to the region of 
t. This gives 



tan a < 



b - b' 



and tan a < 



b' 



3 t 



The induced drag of . one half-wing — is composed ad- 

. 2 

ditively of the drag of half of cone I, W t and the 
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drags of the two half-cones of cone III, and 



r ilIIa (fig. 17). 



nil 



J W i = W Hi + W iIIIx + W HlI 8 ' < 93 > 

Since in cone I, in the plane £ = 0, there is upwash 
velocity, gives a forward thrust which in absolute 

value, however, is smaller than the hack thrust in cone 
III, since the circulation is greater here, We have 

x=t y=b /a 

W ix = P J d x J . c z0 d^y, y x = | - x tan a ) 
x=o y=y a 

where c zo is known from equation (90). We thus have 

2 £=2t/t> d = o 

1 1 - Tl, 8 TT P t 2 K J 

■ ■ e 

In cone I for -1 < 4 < 0: 

11 

1 - Tl 



= / e 2 d I / r U) e ( *) a % (94) 

A l 1 - Tl 1 8 TT p t 2 K 



r= - r 0 K-i 



and therefore 



D 



W iT = 12-( ) / |»d ( / 4 6 (i) d 4 

- - - 1 (1 - lix) 8 - 8- TT • K i l J-. ... -"- J - - - 

£=0 4=_1 

For "briefness we set 

A 

g (*) 1 « - S (i) U = K 3 ( 95a ) 



y g (4) d « = y 1 



k 3 

*=-i. 4=0 



50 - 



E.A.C.A. Technical Memorandum Ko » 897 



-/*g(*)d*-- / « g (*) * K 2 - •- (95b) 

These integrals may be exactly computed. There is obtained 

K 3 = — K a = — (96) 
8 18 



so that finally 



2 



1 (1 - TJj) 2 TT K b 



( ) (97) 



Tha portion ^iju is obtained from equation (94) by 

substituting -g(d) for g($) and taking r = F so 
that 

AAA1 3 1 - 11 r TT K b 

Finally, is obtained by putting in equation (94) 



(i - JCit- ■♦) 



and substituting -g(4) for g(d). By comparison with 
equations (97) and (98) this gives 

W* = W* + IT4 

1 III a a Ii 1 III 1 

and therefore 

\ f t = 3 fi IIXl> + 2 f ijx = 2 W lnii {l + L-} 

^IIIi 
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Substituting .the values from (97) and (98) the induced 
drag of the entire wing is found to be 



4 k 3 pr 0 s 



t tan q ^ 1 3 Kg_ t tan q | 



3 1 - T il tt b . 1-^2 K 3 b 

If, in place of T 0 , there is now substituted the .lift A 
of the entire wing 

a = p b r 0 u 0 



we have 

w ^ _ 16 E3 1 1 / A \ 2 t tan q 

3 tt' (1 - T^Hl + T) x ) 2 P^b u 0 ' . b 

3 K P t tan q 



1 o £.9 u lan q -\ 

ji 1 } 



-2 1 / A \ 3 t tan a 

W. = - 



(1 - TU)(1 + T),) 2 P^O 



(1 - -JU- il- 1^-2 )(ioo) 

L 1 - TJ 8 9.TT b J 



Thus the formula has been found for the induced drag with 
trapezoidal lift distribution. To this must be added the 
wave drag. The latter according to equation (26) and 
table I is 

2(2 + TU) b 1 /JLV (10 l) 

- - - - wave - - X2 - Vu D b/ - U ' 

3(1 + ri x ) s t tan a P 0 

If c w denotes the coefficient of the wave plus induced 
drag then from equations (100) and (101 ) 

— (i ■- — — — X ) (102) 

) 2 ^ 1-T) i 9 TT 7 



/ c a a 4(2+n x ) : 4 X rf f 1 14 

c w / = - + 



4 tan a 3(1+^) 3 (l-*n j) (l+*n J 
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The above f ornrdla differs from the corresponding formulas 
for the rectangular flat plate (equation (80)) and the trape- 
zoidal wing with constant lift distribution (equation (49)) 
in that for small X the induced portion of the drag is 
proportional to X 3 whereas for the other two cases it is 
proportional to X» In figure 20 the coefficient 

/ c a 2 

Cw/ is plotted against the reduced aspect ratio 

f 4 tan a / 

* tan as X for various trapezoid shapes b f /b. It may 
h 

be seen that by far the greatest portion of the drag is 
contributed by the wave resistance. The portion contrib- 
uted by the induced drag, within the range of validity of 
our formulas 9 amounts to a maximum of 11 percent of the 
wave resistance for X =0.5 and b f /b = l/2 # It is 
therefore smaller than for the rectangular flat plate where 
for the same aspect ratio it amounts to 19 percent (fig. 
10). 

VIII. SUMMARY 

UTith the aid of the expressions given by L. Prandtl 
(reference 2) a theory is developed of the airfoil of fi- 
nite span at supersonic speed. As in the case of the 
Prandtl airfoil theory for tho incompressible flow, it is 
a first order approximation theory. Tho airfoil is first 
replaced by a "horseshoe vortex" and the induced velocity 
field of the latter computed. This field is considerably 
different from that of the incompressible flow. From the 
horseshoe vortex there are obtained in the familiar manner 
by superposition more complicated lifting systems. The 
computation of the induced drag, in contrast to the incom- 
pressible case, is for the compressible flow possible only 
if there is first assumed a surface vortex distribution 
and secondly a suitable dropping off of the lift toward 
the wing tips. 

As an example of the "first principal problem" there 
are computed the induced drag and the wing surface shape 
for a wing of trapezoidal plan form with constant surface 
density of the lift. The induced drag, as in the case of 
the incompressible flow, is found to be proportional to 
the square of the lift and depends on the Mach number as 
well as on the aspect ratio. In addition to the frictional 
and induced drag there is present in the supersonic case 
also the wave drag, produced by the sound waves, which 



N.A.C.A. Technical Memorandum Nq. 897 53 



varies as the induced drag. It is therefore only the sum 
of the wave and induced drags that is of practical inter- 
est. 

As an example of the "second principal problem" there 
is computed the lift distribution and induced drag for the 
rectangular flat plate (untwisted rectangular wing). Out- 
side the two Mach cones springing from the leading edges 
of the wing tips the lift density is constant; within 
these cones the lift drops from the full value at the cone 
rim to the value zero at the lateral wing edge. The inte- 
gral equation that arises is independent of the aspect 
ratio and of the Mach number and may he solved numerically 
"by approximate methods. In general for airfoils of normal 
aspect ratios at supersonic flows the greatest portion of 
the total drag is contributed by the wave resistance while 
the induced drag contributes only a small proportional 
part. 

Finally, there is considered the lifting line with 
trapezoidal lift distribution and the lifting surface of 
rectangular plan form whose lift is constant along the 
chord and trapezoidal along the span. For these cases the 
downwash distribution and induced drag are computed. 

Translation by S. Heiss, 
National Advisory Committee 
for Aeronautics. 
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Pigs. 1,8,3,4,5,6 



Figure 3.~ 
Lifting line with 
constant lift 
distribution. 
Downwash distrib- £ 
ution in Mach 
cone. 



Figure 1.- Potential of the 
lifting line. 




Wave of rarefaction^ 



Compression 
shock 




Compression shock v rarefaction 

Figure 4.- Plane sound waves 
at a flat plate. 




Figure 5.- Rectangular 

wing as lifting 
surface with constant lift 
distribution. 



Figure 6.- Rectangular wing as lifting 

surface with constant lift 
distribution. Downwash distribution in 
the wing plane. Continuous curves for 
x< t (at location of wing) dotted curves 
for. x a 2t (behind the wing). 
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Figure 7.- Bectangular wing as Figure 8.- Rectangular wing as 

lifting surface with lifting surface with 



constant lift distribution. constant lift distribution. 

Streamline picture of the y- and s~ Streamline picture of the y- and z- 
velocities in a plane x<t at right Telocity components in the plane 
angles to the axis of the Kach cone. z = 3t at right angles to the axis 

of the toach cone. 




Figure 10.*- Trapetoidal wing with constant lift distribution. 

Coefficients of/ISto wave plus induced drag Cw/ 4tancc 

as a function of ihi> * re due e4 aspect ratio 0 A = t tan<^b for 
various trapezoid shapes 0 »-taiiT/ tanoC. 
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Figure 11.- Trapesoidal wing 

with constant lift 
distribution. Induced downwash 
velocity in section AB (in 
direction of flow) (tanT« 1/V3; 
tanoteV^). 




Pigs. 11,18,13,14 
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Figure 12.- Trapesoidal wing 

with constant lift 
distribution. Induced downwash 
velocity in section CO (at right 
angles to flow direction) 
(tanT«l/V5; tand-V3). 




Figure 14.- Trapesoidal wing 

with constant lift 
distribution. Elevation contour 
lines. (tan<* = V3; tanr « 1/V5). 



Figure 13.- Trapesoidal wing 

with constant lift 
distribution. Profile sections. 
(tanoC = V3; tanT = 1/V5). 
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Pigs, 15,16,17,18 
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Figure 15.- Rectangular plane 
wing. Lift at 

wing edge. 
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Figure 17.- Rectangular surface 

with trapezoidal 
lift distribution. 



Figure 16.- Polars of plane 

rectangular wing 
for various aspect ratios. 




Figure 18,- Lifting line with trapesoidal lift distribution. 

Streamline picture of the y- and s- velocities in 
a plane at right angles to the axis of the Mach cone. 
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Figs. 19.20 




Figure 19.- Rectangular wing as lifting surface with trapesoidal 
lift distributioae Downwash distribution for x<%. 







\1 



its 



4* 



Figure 20.- Lifting surface with trapesoidal lift distribution. 

i c« 

Coefficient of wave plus induced drag Cw/ 4^ tf *• * 
function of the "reduced aspect ratio" A for various values of b'/b. 



